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Abstract We compute the characteristic parameters of magnetic dipole radiation of a 
neutron star undergoing torsional seismic vibrations under the action of Lorentz restor- 
ing force about axis of a dipolar magnetic field experiencing decay. After brief outline of 
general theoretical background of the model of vibration powered neutron star, we present 
numerical estimates of basic vibration and radiation characteristics, such as the oscillation 
frequency, lifetime, luminosity of radiation, and investigate their time dependence upon 
magnetic field decay. The presented analysis suggests that gradual decrease in frequen- 
cies of pulsating high-energy emission detected from a handful of currently monitored 
AXP/SGR-like X-ray sources can be explained as being produced by vibration powered 
magneto-dipole radiation of quaking magnetars. 
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1 INTRODUCTION 

The last two decades have seen an increasing interest in the Soft Gamma Repeaters (SGRs) and 
Anomalous X-ray Pulsars - commonly referred to as magnetars (Duncan & Thompson 1992), the seis- 
mic and radiative activity of which is fairly different from that of rotation powered radio pulsars (e.g. 
Kouveliotou 1999, Harding 1999, Woods & Thompson 2006, Mereghetti 2008, Qiao, Xu & Du 2010). 
The most popular idea is that energy supply of long-periodic pulsating radiation of these fairly young 
neutron stars comes from process involving decay of ultra strong magnetic field. One of such pro- 
cesses could be magneto-mechanical vibrations driven by forces of magnetic-field-dependent stresses 
(Bastrukov et al. 2002). In the development of this line of argument, particular attention has been given 
to torsion vibrations of perfectly conducting stellar matter about magnetic axis of the star under the 
action of magnetic Lorentz force with focus on discrete frequency spectra of toroidal Alfven mode (a- 
mode). In the past, the standing-wave regime of such vibrations has been subject of several investigations 
(Ledoux & Walraven 1958). In works (Bastrukov et al. 2009a, 2009b, 2010) focus was made on the non- 
investigated regime of node-free vibrations in in the static (time-independent) field. The prime purpose 
of these latter works was to get some insight into difference between spectra of discrete frequencies of 
toroidal a-modes in neutron star models having one and the same mass M and radius R, but different 
shapes of constant-in-time poloidal magnetic fields. By use of the Rayleigh energy method, it was found 
that each specific form of spatial configuration of static magnetic field about axis of which the neutron 
star matter undergoes node-free differentially rotational oscillations is uniquely reflected in the discrete 
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frequency spectra by form of dependence of frequency upon overtone £ of vibrations. The subject of 
present our study is radiative activity of quaking neutron star powered by energy of Alfven vibrations 
in its own magnetic field experiencing decay. Part of this project has been reported in recent workshops 
and conferences and in short paper (Bastrukov et al. 201 1). Taking this into account, in this article only 
a brief overview of theoretical background of the model is given. The focus is placed on numerical com- 
putation of basic characteristics of vibration and radiation which are of interest in observational search 
for such objects. 



2 GENERAL BACKGROUND OF THE VffiRATION POWERED NEUTRON STAR MODEL 

Only a brief outline of this model is given here and more details can be found elsewhere (Bastrukov 
et al. 201 1). The Lorentz-force-driven differentially rotational vibrations of perfectly conducting matter 
of neutron star about axis of poloidal internal and dipolar external magnetic field evolving in time are 
properly described in terms of material displacements u obeying equation of magneto-soUd-mechanics 

p{r) u(r, = -L [V X [V X [u(r, t) x B(r, t)]]] x B(r, t) (1) 
in 

u{r,t) = [uj{r,t)xr], a;(r, t) = ^t[VxW] (2) 

The field u(r, t) is identical to that for torsion node-free vibrations restored by Hooke's force of elastic 
stresses (Bastrukov et al. 2007, 2010) with x(r) = Ag fi{r) Pi{cofi 9) where fi{r) is the nodeless 
function of distance from the star center and Pe{cos9) is Legendre polynomial of degree £ specifying 
the overtone of toroidal mode. In the amplitude a{t) is the basic dynamical variable describing time 
evolution of vibrations which is different for each specific overtone. The bulk density can be represented 
in the form p{r) = p<P{r), where p is the density at the star center and (f){r) describes the radial profile 
of density which can be taken from computations of neutron star structure relying on realistic equations 
of state accounting for non- uniform mass distribution in the star interior (e.g., Weber 1999). The central 
to the subject of our study is the following representation of the time-evolving internal magnetic field 
B(r,i) = B{t)h{r), where B{t) is the time-dependent intensity and b(r) is dimensionless vector- 
function of the field distribution over the star volume. The gist of the energy variational method of 
computing frequency of nod-free Alfven vibrations consists in the following separable representation of 
material displacements u(r, t) = a(r) a{t). Scalar product of (??) with this form of u and integration 
over the star volume leads to equation for amplitude a{t) having the form of equation of oscillator with 
time-depended spring constant 

Ma{t) + !Cit)a{t) = 0, (3) 
A4 = pm, m = 



IC — — k, k 
47r 



J <^(r) a(r) • a(r) dV, a = At\7 x [r Mr) Pe {cos 9)] 
J a(r) • [b(r) x [V x [V x [a(r) x b(r)]]]] dV 



As for the general asteroseismology of compact objects is concerned, the above equations seems to be 
appropriate not only for neutron stars but also white dwarfs (Molodtsova et al 2010, Bastrukov et all 
2010a) and quark stars. The superdense material of strange quark stars is too expected to be in solid 
state (Xu 2003, 2009). In Fig.l we plot the fundamental frequency va = WA/27r (where uja = va/R) 
and the period Pa = '^a^ of global Alfven oscillations 
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as functions of intensity of constant in time magnetic field B for solid star models with masses and radii 
of typical white dwarfs, neutron stars and quark stars. 




Fig. 1 The basic frequencies and periods of global Alfven oscillations as functions of mag- 
netic field for the typical mass and radius of white dwarfs, neutron stars and quark stars. The 
difference in mass and radius can be regarded as reflecting difference in underlying equations 
of state of matter in these compact objects. 



The total vibration energy and frequency are given by 



EaH) ^ ^ ^ — Lu(t) ^ \ ^ B(t)K, s. (5) 

Here AI and R are the neutron star mass and radius and s is the parameter depending on overtone of 
Alfven toroidal mode and the depth of seismogenic layer. The energy conversion of above magneto- 
mechanical vibrations into magneto- dipole radiation which is governed by equation 

^^-T'W, Pit)^^^Si.'{t). (6) 

The axisymmetric torsional oscillations of matter around magnetic axis of the star are accompanied by 
fluctuations of total magnetic moment preserving its initial (in seismically quiescent state) direction: 
fi = fin = constant. The frequency a;(t) of such oscillations must be the same for both fluctuations of 
magnetic momenta 6fj,{t) and above magneto-mechanical oscillations which are described in terms of 
a{t), namely 

Sil{t)+uj^{t)5fj,{t) ^0, a(t) +w2(t)a(t) = 0, uj^ {t) ^ (t) . (7) 

This suggests 6fi{t) = fia{t). On account of this the equation of energy conversion is reduced to the 
following law of magnetic field decay 

dB(t) , 2n^K^ 

-^--iB^t), ,= ^ . constant, (8) 
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Fig. 2 The figure illustrates the effect of magnetic field decay on the vibration frequency and 
amplitude of quadrupole toroidal a-mode presented as functions of x ^ t/r. 

S(i) = -£E=, r-i=27i?^(0). (9) 

vi + 

Thus, the magnetic field decay resulting in the loss of total energy of Alfven vibrations of the star causes 
its vibration period P ^ B^^ to lengthen at a rate proportional to the rate of magnetic field decay. 
The time evolution of vibration amplitude (the solution of equation for a{t), obeying two conditions 
a{t = 0) = ao and a{t = r) = 0) read^](e.g., Polyanin & Zaitzev 2004) 

a{t) = Cs^'^{,h{z{t))-r]Yi{z{t))}, z{t)^2l3s^l^{t), s = 1 + t/r (10) 

where Ji(z) and Yi(z) are Bessel functions (e.g., Abramowitz & Stegun 1972) and 

2 2Ea{Q) 2Ea{Q) , R'iO) 

Here by Ea{0) is understood the average energy stored in torsional Alfven vibrations of magnetar. If 
all the detected energy -Eburst of ^-ray outburst goes in the quake-induced vibrations, E'burst ~ E^., 
then the initial amplitude uq is determined unambiguously. Thus, the magnetic field decay is crucial to 
the energy conversion from Alfven vibrations to electromagnetic radiation whose most striking feature 
is the gradual decrease of frequency of pulses (period lengthening). The magnetic-field-decay induced 
lengthening of period of pulsating radiation (equal to period of vibrations) is described by 

p«)^p(o)|l+«Ml■'^ pw-^ n^;';4.,. , 

p2(o) ^ ^ 27r 

^(0) = (12) 



2P{t)P{ty ' ' kB{0) 



The authors are indebted to Alexey Kudiyashov (Saratov, Russia), for the helpful assistance in solving this equation. 



Vibration powered neutron star 
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Fig. 3 Artist view of quadrupole overtone of torsional Alfven vibrations of magnetar about 
magnetic axis producing oscillations of lines of dipolar magnetic field defining the beam 
direction of outburst X-ray emission. 

On comparing r given by equations (|9]l and ( fT2] i. one finds that interrelation between equilibrium equi- 
librium value of the total magnetic moment of a neutron star of mass M = 1.2 Mq and radius i? = 15 
km vibrating in quadrupole overtone of toroidal a-mode, pictured in Fig. 3, is given by 

/I = 5.5 X 10^^ P[t)P{t), Gcm^. (13) 

This illustrative picture shows that pulsating radiation in question belongs to the periodic changes 
of polarization of magneto-dipole radiation powered by Alfven seismic vibrations of neutron star. In the 
reminder we present numerical estimates of characteristic parameters of of pulsating magneto-dipole ra- 
diation produced by quaking neutron star at the expense of energy of Lorentz-force driven differentially 
rotational vibrations about axis of dipole magnetic moment. 

3 NUMERICAL ANALYSIS 

To get an idea of the magnitude of characteristic parameter of vibrations providing energy supply of 
magneto-dipole radiation of a neutron star, in Table 1 we present results of numerical computations of 
the fundamental frequency of neutron star oscillations in quadrupole toroidal a-mode and time of decay 
of magnetic field r as functions of increasing magnetic field. 

As was emphasized, the most striking feature of considered model of vibration powered radiation 
is the lengthening of periods of pulsating emission caused by decay of internal magnetic field. This 
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Table 1 The Alfven frequency of Lorentz-force-driven torsion vibrations, i^a, and their life- 
time equal to decay time of magnetic field, t, in neutron stars with magnetic fields typical for 
pulsars and magnetars. 





M(Mq) 


K(km} 




!^a(Hz) 


r(yr) 


Pulsars 


0.8 


20 




3.25 X 10^'' 


4.53 X 10'" 




1.0 


15 




2.52 X 10"^ 


2.98 X lO'^ 


Magnetars 


1.1 


13 


10^* 


0.22 


7.4 X 10^ 




1.2 


12 




2.06 


1.31 




1.3 


11 




1.89 


2.38 




1.4 


10 


10'' 


17.4 


4.44 X 10"" 



suggests that this model is relevant to electromagnetic activity of magnetars - neutron stars endowed 
with magnetic field of extremely high intensity the radiative activity of which is ultimately related to 
the magnetic field decay. Such a view is substantiated by estimates of Alfven frequency presented in 
the table. For magnetic fields of typical rotation powered radio pulsars, B ^ 10^^ G, the computed fre- 
quency VA is much smaller than the detected frequency of pulses whose origin is attributed to lighthouse 
effect. In the meantime, for neutron stars with magnetic fields B ^ 10^^ G the estimates of va are in 
the realm of observed frequencies of high-energy pulsating emission of soft gamma repeaters (SGRs), 
anomalous X-ray pulsars (AXPs) and sources exhibiting similar features. According to common belief, 
these are magnetars - highly magnetized neutron stars whose radiative activity is related with magnetic 
field decay (e.g.. Woods & Thompson 2006). The amplitude of vibrations is estimated as 

1/2 

= 3.423 X 10^'' e]/IqB^IR^^^\ (14) 

where Eaao ~ Ea/{1Q*'^ erg) is the energy stored in the vibrations. Re = i?/(10^ cm) and B14 = 
5/(10^'* G). The presented computations show that the decay time of magnetic field (equal to duration 
time of vibration powered radiation in question) strongly depends on the intensity of initial magnetic 
field of the star: the larger magnetic field B the shorter time of radiation t at the expense of energy of 
vibration in decay during this time magnetic field. The effect of equation of state of neutron star matter 
(which is most strongly manifested in different values mass and radius of the star) on frequency va is 
demonstrated by numerical vales of this quantity for magnetars with one and the same value of magnetic 
field B = 10^^ G but different values of mass and radius. 

According to above presented analytic results, the neutron star luminosity powered by neutron star 
vibrations in quadrupole toroidal a- is given by 

„2 

V ^K-^B-^{t)a^{t), n = (l/2)B(0)i?^ B{t) ^ B{0)[1 + t/r]-^/^, (15) 

a{t)=Cs^/^{Ji{z{t))-rjYi{z{t))}, z ^ 2lj{0)t (1 + t/r) (16) 

The presented in Fig. 4 computations of power of magneto-dipole radiation of a neutron stars (of one 
and the same mass and radius but different values of magnetic fields) exhibit oscillating character of 
luminosity. The frequency of these oscillations equal to that of torsional Alfven seismic vibrations of 
neutron star The practical usefulness of presented computations is that they can be used as a guide in 
observational quest of vibration powered neutron stars among the currently monitoring AXP/SGR-like 
sources. 

4 SUMMARY 

It is generally realized today that the standard model of inclined rotator, lying at the base of our un- 
derstanding of radio pulsars, faces serious difficulties in explaining the long-periodic (2 < P < 12 s) 
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Fig. 4 Time evolution of luminosity of magneto-dipole radiation powered by energy of tor- 
sional Alfven seismic vibrations of neutron star with mass M = 1.2Mq and radius i? = 15 
km with intensities of magnetic field. 



pulsed radiation of soft gamma repeaters (SGRs) and anamalous X-ray pulsars (AXPs). Observations 
show that persistent X-ray luminosity of these sources (10"^'* < Lx < 10^^ erg s~^) is appreciably 
(10-100 times) larger than expected from neutron star deriving radiation power from energy of rota- 
tion with frequency of detected pulses. It is believed that this discrepancy can be resolved assuming 
that AXP/SGR-like sources are magnetars - young, isolated and seismically active neutron stars whose 
energy supply of pulsating high-energy radiation comes not from rotation (as is the case of radio pul- 
sars) but from different process involving decay of ultra strong magnetic field, lO^'' < B < 10^® G. 
Adhering to this attitude we have investigated the model of neutron star deriving radiation power from 
the energy of torsional Lorentz-force-driven oscillations in its own magnetic field experiencing decay. 
It worth noting that such an idea is not new and first has been discussed by Hoyle et al. (1964), still 
before the discovery of pulsars (e.g., Pacini 2008). What is newly disclosed here is that such radiation is 
possible when magnetic field is decayed. Since the magnetic field decay is one of the most conspicuous 
features distinguishing magnetars from rotation powered pulsars, it seems meaningful to expect that at 
least some of AXP/SGR - like sources can be vibration powered magnetars. Working from this, we have 
presented extended numerical analysis of the model of vibration powered neutron star whose results 
can be efficiently, as is hoped, utilized as a guide for discrimination of vibration powered from rotation 
powered neutron stars. The application of presented here analytic and numerical results to analysis of 
specific currently monitoring sources will be be presented in forthcoming papers. 
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